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ABSTRACT 

We present numerical calculations of the charge and statistics, as extracted from Berry 
phases, of the Laughlin quasi-particles, near filling fraction 1/3, and for system sizes of 
up to 200 electrons. For the quasi-holes our results confirm that the charge and statistics 
parameter are e/3 and 1/3, respectively. For the quasi-electron charge we find a slow 
convergence towards the expected value of —e/3, with a finite size correction for N electrons 
of approximately — 0.13e/A^. The statistics parameter for the quasi-electrons has no well 
defined value even for 200 electrons, but might possibly converge to 1/3. Most noteworthy, 
it takes on the same sign as for the quasi-holes, due to terms that have previously been 
ignored. The anyon model works well for the quasi-holes, but requires singular two-anyon 
wave functions for modelling two Laughlin quasi-electrons. 
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1 Introduction 



It is now widely accepted that the fractional quantum Hall effect arises due to quasi-particle 
excitations created when the filling fraction of the lowest Landau level moves away from 
its preferred values. The quasi-particles bind to impurities and thereby ensure that the 
plateaus in the conductivity are formed. Laughlin examined the filling fractions 1/m, with 
m an odd integer, and argued that the quasi-particles have charge ±e/m (— e is the electron 
charge) He showed that these values for the charge imply a conductivity plateau at 
the value e^/(m/i), as observed. He also offered explicit trial wave functions describing the 
ground state as well as the quasi-particle excitations for the 1/m case. Later, Haldane p| 
and Halperin |^ examined the hierarchy structure of the Hall states, and suggested that 
the quasi-particles obey fractional statistics, i. e. that they are anyons . 

Arovas, Schrieffer and Wilczek derived the charge and statistics of the Laughlin quasi- 
holes in Ref. They examined the Berry phase corresponding to one quasi-hole 
encircling the origin and interpreted this as an Aharonov-Bohm phase 0. The charge 
was then found to be e/m, in confirmation of Laughlin's result. They also considered a 
pair of quasi-holes encircling one another, and related the two-particle contribution of the 
Berry phase to the quasi-hole charge. Interpreting this two-particle contribution as an 
anyon interchange phase they found the anyon statistics parameter to have the value 1/m, 
equal to the fraction of the elementary charge. The Laughlin quasi-electrons were examined 
along the same lines in Ref. 0. Within the approximations used, the results imply that 
the charge and statistics parameter have the values —e/m and —1/m, respectively. 

The statistics satisfied by the quasi-particles in the quantum Hall system was also 
examined in Ref. |^, where the exclusion statistics parameter |jlO| was considered. The 
results did not rely on any specific trial wave function but rather on state counting based 
on numerical simulations for interacting electrons on a sphere. The value of the one- 
dimensional exclusion statistics parameter was found to be 1/m in the case of quasi-holes, 
while the value 2 — 1/m was found for the quasi-electrons. Both cases were examined near 
to the magic filling fraction 1/m with m = 3. 

The exclusion statistics parameter is in principle the same parameter as one reads 
from the Berry phase, although with an opposite sign for the quasi-electrons, because 
their charge is negative. The relation is briefly discussed in the Appendix below. Thus 
the predictions for the anyon parameter, based on the numerical results for the exclusion 
statistics parameter, would be 1/m for quasi-holes and —2 + 1/m for quasi-electrons. The 
values 1/m and —2 + 1/m for the anyon parameter define of course the same particle 
statistics, but we distinguish between them here in the way we define the correspondence 
between quasi-particles and anyons. Thus, in the case of quasi-electrons, 1/m and — 2-|-l/m 
would represent the same species of anyons, but different anyon states, the 1/m state 
singular and the —2 + 1/m state non-singular. 

It is interesting to note that the numerical results for the exclusion statistics parameter 
of realistic quasi-holes and quasi-electrons are easily interpreted in terms of anyons of 
positive and negative charge, respectively, with non-singular wave functions, and with the 
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same statistics in the two cases. That quasi-holes and quasi-electrons should have the same 
statistics, is also what one would expect if one regards them as antiparticles of each other. 
The anyon representation for the quasi-particles, first suggested by Laughlin [TTI, W2 



was considered in more detail in Ref. [|T3|. A careful examination of the quasi-electron case 
revealed that the Berry connection calculated from the Laughlin wave function actually 
contains some terms not earlier considered in the literature. These extra terms arise be- 
cause the inverse of the operator used to create a quasi-electron is not simply the complex 
conjugate of the inverse quasi-hole operator. There is no obvious argument for neglecting 
them, and this motivates a closer examination of whether these wave functions really rep- 
resent excitations carrying the charge and statistics parameter characterizing the physical 
quasi-electrons in the quantum Hall system. There are no similar problems associated with 
the Laughlin quasi-hole wave functions. 

The purpose of the work presented here was to compute the Berry phases for the 
Laughlin quasi-particles in order to derive their charge and statistics, following the ideas of 
Arovas et al. . The calculations were performed both for quasi-holes and quasi-electrons, 
with special attention to the latter, at the filling fraction 1/3. For a small number of 
electrons (< 8) we could integrate exactly, but for larger and more realistic systems we 
had to use Monte Carlo integration with importance sampling according to the Metropolis 



algorithm [11 



In the case of a single quasi-hole as well as for a pair of quasi-holes, we find that the 
Berry phases computed for systems with up to 200 electrons behave as expected. Thus the 
charge e/3 and the statistics parameter 1/3 for the quasi-holes are both confirmed. 

For a single quasi-electron the extra terms affect the Berry phase so as to imply a charge 
that for a system with up to 200 electrons is slightly below —e/3, and also the boundary 
effects are larger than for the quasi-hole. For a pair of quasi-electrons, the computations 
give a definite short distance behaviour of the statistics phase, but at distances larger than 
about one magnetic length, which is roughly the size of a quasi-electron, the phase does 
not settle down to an asymptotic value, as one would like it to. Moreover, to the extent 
that the statistics parameter is a meaningful quantity, it has the same sign as for the quasi- 
holes, in clear contradiction to earlier results for the Laughlin quasi-particles. This sign 
would correspond to singular two-anyon wave functions, and hence a negative exclusion 
statistics parameter. 

Within the last year there has been reported direct experimental evidence for the frac- 



tional charge of the quasi-particles in the fractional quantum Hall system |]14[. According to 



a conjecture proposed in Ref. [T^, fractional charge implies fractional statistics. Thus our 
investigation of the charge and statistics of the Laughlin quasi-particles has some current 
interest for comparison with the real quantum Hall system. 
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2 Charge and Statistics From Berry Phases 



Consider a system of particles (electrons) of charge — e (e > 0), moving in the {x,y) 
plane, with a magnetic field B = Bbz {B > 0), perpendicular to the plane. The position 
of a particle may be described by the complex coordinate z = (x + iy)/{^/2l), where 
/ = 1/^/eB is the magnetic length. We set h,c= 1. 

Due to the magnetic field there is an Aharonov-Bohm phase associated with the prop- 
agation of a charged particle around a closed loop [0. If the path is a circle of radius 
p = ^/2lr, right handed relative to the direction of the external magnetic field, and if the 
particle has charge q, this phase is given by 

7 = vrp^fig = 27r^r2 . (1) 

Consider a quasi-particle (quasi-electron or quasi-hole) excitation, of charge q, localized 
at the position zq and described by the normalized wave function 

Zn) = {Zi,..., Zn I Zq) . (2) 

If this quasi-particle is moved around a closed loop, there arises the so called Berry phase, 
which is the integral along the path of the Berry connection 0. Let the path be a circle 
around the origin, parametrized as zq = re'*^ with 4> running from to 27i. The Berry 
connection is then defined by 

^ = i (zol <9<^ ko) • (3) 

The charge q of the quasi-particle is now determined by setting the Berry phase equal to 
the Aharonov-Bohm phase corresponding to the same path In this way q is expressed 
in terms of parameters of the fundamental particles (electrons) in the system. In a finite 
system the charge defined in this way will depend on the distance r from the origin, but 
it is reasonably well defined as long as the r dependence is small for circles that are well 
inside the system. 

Suppose the quasi-particle is described by a normalized iV-particle state of the form 

-| oo 

l^o) = ^EVaHO- (4) 

V-'l z=o 

Here ai are expansion coefficients, |/) are orthonormal basis states, and Ji is introduced for 
normalization, 

oo 

h = T.r'W- (5) 

/=o 

Then the Berry connection depends on r but not on (p, so that the Berry phase is /5i(27r) = 
2tt dj3i/d(j). The expression for the charge becomes especially simple, and is given by the r 
dependence of Ji as 

q 1 dpi d 



e d(t> dr'' 



In/i. (6) 
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Now suppose there are two quasi-particle excitations simultaneously, located symmet- 
rically about the origin at the positions ±zo- The parametrization zq = re"^ then describes 
a counterclockwise interchange of the quasi-particles if we let run from to vr. If we 
assume that the two quasi-particle state is described by a state analogous to (^), 

oo 

ko,-^o) = ^E-^o'fe/IO , (7) 

V-f 5 ;_r, 



'2 1=0 

then a Berry connection dj32/d(j) corresponding to this interchange may be defined in the 
same way as (H), 

= Hzo^—Zo\d^\zo,—zo) . (8) 
We subtract the single-particle contributions due to the magnetic field, and define 

-iy=- - 2A(7r)) = ^{f32- 2A) = -r^ A (inJ^ - 21nJ0 . (9) 

TT d0 dr^ 

We will refer to u as the "anyon parameter", since the residual Berry phase — i/tt can be 
identified with minus the statistics phase of the particles. 

As a technical point, note that Eq. gives u as the small difference of two quantities 
that are large for large r. When computing u by the Monte Carlo method, it is important 
to evaluate both integrals Ji and I2 simultaneously, using the same sample of random 
points, since the statistical errors will then be correlated and tend to cancel. 

The Laughlin quasi-particles of the quantum Hall system are described by states of the 
same form as in Eq. (^, or by states of this form with zq replaced by Zq- In the next section 
we give the explicit expressions for the wave functions. For now it suffices to know that 
the expansion states corresponding to |/) in (D are orthogonal, and the same is true for 
|/') in d^), which are different from |/). A non-normalized quasi-electron state is expanded 
as a polynomial in zq, whereas a quasi- hole state is expanded in Zq instead of zo- The 
complex conjugation shows up as a difference in sign in the relation between the Berry 
connection and the normalization factors Ji and l2- We use the notation I^^^, J^q^jiqe^ j2qe 
for the normalization factors for one and two quasi-holes and one and two quasi-electrons, 
respectively, with N electrons in the system. The following expressions for the charges and 
anyon parameters of these quasi-particles then result. 



e dr^ 

^ ln/i<^% (11) 



e dr^ 

d 



= -r'— (ln/2^^-21n/^^^) , (12) 
dr"' 



^qe _ r^— (ln/2q"-21n/i^n . (13) 
dr^ 
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This shows immediately that the charge is positive for the quasi-holes and negative 
for the quasi-electrons. The sign of the anyon parameter in the two cases is not obvious 
directly from these formulae. As discussed in the Appendix, the anyon model predicts that 
> and z/^*^ < 0, when only non-singular two-anyon wave functions are considered, 
whereas the ranges — 1 < u'^^ < and < u^'^ < 1 correspond to singular two-anyon wave 
functions. 



3 Computational methods 



One quasi-hole 

The non-normalized A^-electron wave function for one quasi-hole at the position Zq is 

N N 

^l^\z, z*) = i,. A*™ n(^r - ^o) = ^0 A*'" Y.(-zlfSuA^*) . (14) 

i=l k=0 

We simplify our notation by writing z = {zi,...,zn). Here 1/m is the filling fraction, ipo is 
the harmonic oscillator ground state wave function, 

^0 = e-^E,=il^«l% (15) 

and A is the Vandermonde determinant, 

A = ll{z,-Zk). (16) 

j<k 



The elementary symmetric polynomials Sk,N [0 are implicitly defined by the above ex- 



pansion, which brings the quasi-hole state explicitly to the form shown in Eq. (|^), except 
for the substitution zq ^ Zq. For computations we use the following recursion relation, 

Sk,N{z) = Zn Sk,N-l{z) + Sk-l,N-l{z) , (17) 

with 5*0,0 = 1 and all other 5*^,0 = 0. 

The normalization integral, which determines the charge according to Eq. ( TH ) , is given 

by 

Ji^^r^) = Id^^^lvj/J^p (18) 

= fd'^^z^o'lArfilz.^zol' (19) 
■' k=i 

= Y.r''' [d'''z^p,'\Ar\S,M^)\' (20) 

k=0 

N 

= Er^'ll'"'. (21) 

k=0 
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Since we want the logarithmic derivative with respect to r^, we need compute only the 
ratios We are able to do this analytically only for very small N. For larger N 



we do Monte Carlo integration by the Metropolis algorithm |rT[. This method works well 
in particular because 

l]f = j d^^zi;o^\A\^^"' , (22) 

and the integration measure (^^z iJjq^ | Ap™ is common to all the integrals l]^- We interpret 
V'o^ [Ap"* as a Boltzmann factor e~^^ of classical statistical mechanics, thinking of (3 as 
the "inverse temperature" and V as the "potential energy" , with 

(3V = Y.\z^\^ -2mY^\Ti\zi- Zj\. (23) 

i i<j 

The classical system is a plasma with two-dimensional Coulomb repulsion, in an external 
harmonic oscillator potential |]1T|. The factor H^i \zk — -ZoP, which we do not include 



in the integration measure, represents in this plasma analogy a Coulomb repulsion of the 
electrons from the hole position zq. 

We compute the ratio 1^°"^ /l]^^ as the time average of \Sk,N{,z)\'^ over a time sequence of 
A^-electron configurations generated by the following random dynamics. We loop through 
all electrons, and for the i-ih electron we generate a random "trial jump" Zi ^ Zi + Azi in 
such a way that /\zi and —/S.Zi are equally probable. The trial step is accepted or rejected 
depending on whether the Boltzmann factor e~^^^, where AV^ is the change in potential 
energy, is larger or smaller than a random number generated uniformly between and 
1. This procedure gives the desired distribution of configurations, because the dynamics 
satisfies the principle of detailed balance: the ratio of probabilities for jumping from A to 
B or from i? to A is the Boltzmann factor e"^*^^^'^^"^*^'^-*-'. 

One quasi-electron 

Now consider a quasi-electron at the position zq. Laughlin's proposed wave function is a 
polynomial in zq, 

/ N \ N 

v^lf (^, z*) = V^o n(^< - ^o) A*- = ^0 E(-^o)'5fc,^(5*) A*™ , (24) 

\j=l / fc=0 

where d* = {dz*,...,dz*^). The normalization integral can be written in several different 
ways. We may expand it as a polynomial in r^, 

Ji^<=(r2) = J d''^z\^l'^'\^ (25) 

TV 

J2r'' Jd'''ziJo'\SkAd)An' (26) 



k=0 

N 

„2k 7-lqe 



Y.r''ir, (27) 

fc=0 
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or rewrite it by partial integration as 



N 



Jlqe(^2) ^ f ^2N^^^2 \Ari[{\z, - Zo\' - I) . (28) 

k=l 

The last form shows clearly the difference between the quasi-hole and quasi-electron cases. 
If we could neglect the —1 in each factor — zop — 1, then the integrals would be identical, 
and the charges would be the same, just with opposite signs. This approximation seems 
hard to justify, since the average number of electrons within unit distance from an arbitrary 
point Zo, assuming constant density within the electron droplet, is vr/m, which is close to 
one if m = 3. A more likely possibility would be that the normalization integral does 
indeed change, but only by a factor which is mainly independent of r, and which therefore 
does not change the Berry phase. 

Note that the integrand in Eq. (^) may even be negative, in spite of the fact that the 
integrand of the original normalization integral is explicitly non-negative. This can happen 
because the derivation of Eq. ( pS]) involves 2N partial integrations. The cancellation be- 
tween positive and negative contributions will cause some loss of precision when this form 
of the integral is used for numerical evaluation. 

Before proceeding we want to comment on the origin of the extra —1. Remember that 
the wave function of one particle in a magnetic field, in the lowest Landau level, has the 
form f{z*) e~2l^l , where / is an analytic function. The projection onto the lowest Landau 
level of the operator z is the operator d^* acting on the space of functions analytic in z* 



Tq . This is seen by partial integration, which gives that 

d'ze-^''\g{z*)rd.,f{z*) = J d^ze'^^* {g{z*)yzf{z*) (29) 

for two analytic functions / and g. Partial integration also gives that 

I d'ze-^^'iid,* - z,)g{z*)nd,, - Zo)fiz*) = J d^ze— *((7(^*))* i\z - Zo\' - l)/(^*) .(30) 

In the projection onto the lowest Landau level there is an operator ordering problem, 
because of the fundamental fact that [z*,(92*] = —1. Thus, the factor |2; — 2;op in the 
quasi-hole normalization integral corresponds to the ordering {z — zq){z* — Zq), whereas the 
factor \z — ZqI"^ — 1 in the quasi-electron integral corresponds to {z* — Zq){z — zq). 

We may use Eq. (|28|) to compute the integrals ll^^ defined in Eq. (|26|) , by expanding 

N N 

Y[{\z^ - Zo\' - I) = T.<i^o'z*^. (31) 

n=l k,l=0 

The coefficients a^i satisfy the following recursion relation, 

with ttg Q = 1 and all other a° ^ = 0. Only the diagonal coefficients give nonzero 
contributions to the integral, but in order to compute them we have to compute some, 
though not all, of the off-diagonal coefficients Note that is real, and more generally, 
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Two quasi-holes 



For two quasi-holes at the positions Za and zi, the wave function is 

N 

= ^0 A*- Wiz* - zDiz* - zl) . (33) 



i=l 



The normahzation integral also for this state takes the form of a probability distribution 
for a system with two-dimensional Coulomb interactions, 

N 

lit. = / d'^'^^o' n - ^a){zk - , (34) 

•' k=l 

which is analogous to Eq. (p!9D. If we let the quasi-holes be located symmetrically around 
the origin, Za = —Zb = zq, then the state is on the form of Eq. (^, with zq ^ z^, and the 
normalization integral can be expressed as 

N 
A:=0 

Here 

J^^ = / d'^'z^o' \S,A4,-,4)\' (36) 



are again normalization integrals for certain angular momentum eigenstates. The elemen- 
tary symmetric polynomials now depend on z"^ instead of z. Numerical calculation of these 
integrals is just as straightforward as in the case of a single quasi-hole, since the recursion 
formula in Eq. ([T7|) is valid with the substitution z ^ z^. 



Two quasi-electrons 

Two quasi-electrons at Za and Zb are described by the wave function 

= ^0 (n (5.. - z^) (9.. - Zb)^ A™ , (37) 
which yields, by partial integration, the normalization integral 

llZ, = / d'"^^ V-o' |A|2™ n (kfc - Za?\Zk - Zb\^ - \2Zk -Za- Zbf + 2) . (38) 

•' k=l 

Comparing with the quasi-hole integral in Eq. (|3^, we see that for two quasi-electrons 
there are correction terms, like in the case of a single quasi-electron. 

We choose to calculate Berry phases for two quasi-electrons located symmetrically 
around the origin, with Za = —Zb = zq, which again ensures that our state is of the 
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form given in Eq. (^. Tlie normalization integral needed to find the statistics parameter 
according to Eq. (|13D , takes the form 



N 

E 

A:=0 



J2qe(^2)^^^4.^2cie^ (39) 

where 



= / d^^.^o^ \s.Adl,-,dl)A-\^ . (40) 

However, it is again simpler to compute ll^^ by first integrating partially to eliminate the 
derivatives and then expanding 

n(k„^-VP-4|.„P + 2)= 5: b^,Zo''zf. (41) 

n=l k,l=0 

The coefficients b^i satisfy the recursion relation 

bl, = {\zX - 4\z^\' + 2) b-f - <2 b-zl - z^' b-J^, + b-zl_^ , (42) 

with 6qq = 1 and all other 6°; = 0. Again it is only the diagonal coefficients b^j. that give 
nonzero contributions to the integral. 



4 Results 

In this section we will present results from the Monte Carlo calculations described above. 
We have plotted, as functions of the dimensionless radius r of a circular loop, the charge 
and statistics parameter for Laughlin quasi-holes and quasi-electrons, extracted from Berry 
phases according to Eqs. (0), (0) and (|T3|). 

The case we have studied is m = 3, i.e. 1/3 filling of the lowest Landau level. For the 
quasi-holes we find Berry phases that imply the charge e/3 and the statistics parameter 
1/3, in confirmation of the results from Ref. [^. 

For the quasi-electrons our results are not nearly as unambiguous. The difference as 
compared to the quasi-hole case is due to the extra terms in the normalization integrals, 
as discussed above, and for the first time in Ref. |jl3|. These terms affect the value of the 
quasi-electron charge, making it slightly larger than e/3 in absolute value. Nevertheless 
the charge is close to the expected value of —e/3, and it seems plausible that the deviation 
is a finite size effect, vanishing in the limit of infinitely many electrons. 

The statistics parameter of the quasi-electrons is where the difference shows up most 
clearly. In particular, the statistics parameter is found to have the same sign as for quasi- 
holes, and not the opposite sign as would be the case if the extra terms were unimportant. 
Furthermore, as seen in the figures presented below, there is no clear indication in our 
numerical results that the statistics parameter will converge either to 1/3 or to any other 
value as the number of electrons increases. 
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Consider first the bulk value of the quasi-hole charge q'^^ as given by Eq. (]I0|). In Fig. |1| 
the quantity q'^^/e is shown as a function of r, the dimensionless distance from the origin 
to the quasi-hole. For filling fraction 1/m the radius of a circular electron droplet with N 
electrons is approximately ymN . With m = 3 and for the values of considered here 
this means that the droplet boundaries are at r 7.7, 12.2, 15.0, 17.3 and 24.5. Fig. |l| 
then clearly illustrates that as long as the quasi-hole is well inside the droplet, its charge 
is well defined and equal to 1/3 of the absolute value of the electron charge. This result 
is in accordance with the bulk value proposed in Ref. and later calculated in Ref. 0. 
We also note that the boundary effects weaken as the number of electrons grows. 

In Fig. ^ we show a more detailed picture of the 100 electron case, where we have 
collected a statistics of 409 million points in the Monte Carlo integration. The bulk value, 
which we define as the value for not too small and not too large radius, is seen to take the 
expected value of 1/3 to at least four decimal digits. We consider the deviations from 1/3 
to be not statistically significant, this we checked in some cases by splitting the data in 
two or more parts and then plotting curves for each partial data set. The deviations from 
1/3 at large r are clearly edge effects due to the finite system size. 

The deviations at small r we believe to be also finite size effects that can be understood. 
In fact, the logarithmic derivative of /^^'^(r^) is determined at r = by only two coefficients 
1^°^^ and l]^, whereas for medium values of r it is a kind of average over many coefficients. 
Each coefficient may be subject to fluctuations, due to the finite Monte Carlo statistics 
and perhaps also due to the finite system size. Such fluctuations are more hkely to be 
averaged out in quantities depending on many coefficients. 

Fig. P shows the charge of the quasi-electron in terms of the absolute value of the 
electron charge, g'^'^/e, which is given by Eq. ([TT|) and is expected to have the value —1/3. 
As for the quasi-holes, the charge is well defined in the bulk, and the boundary effects 
decrease with increasing electron number. However, the bulk value is slightly lower than 
—1/3. This result fits together with the common belief expressed elsewhere in the literature 
only if the deviation vanishes with increasing electron number, which may well happen, 
although our present data are not conclusive evidence. Comparing Figs. ^ and ^ we see 
that the edge effect due to the finite size of the electron droplet is larger in the case of 
a quasi-electron than it is for the quasi-hole. Both of these features, the deviation from 
—1/3 and the larger edge effects, are due to the extra terms in the normalization integrals. 

Fig. ^ shows q"^^ /e for small distances, r < 6, including also curves computed for 100 
and 200 electrons. The results shown may be parametrized roughly as 

^ = (43) 

e 3 ^ ' 

To the extent that this parametrization is valid, the conclusion is that the deviation from 
— 1/3 is a finite size effect. 

We focus next on the statistics parameter for the Laughlin quasi-holes, Eq. ([T^). 
Fig. ^ shows this quantity as a function of r, half the distance between the two holes. When 
the two quasi-holes are well inside the electron droplet and also reasonably well separated. 
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the quantity u'^^ is seen to have a rather well defined constant value of 1/3. These two 
characteristics together make it meaningful to state that the quasi-holes have the anyon 
statistics parameter 1/3. Another general feature is that when the two quasi-holes are 
close to the edge of the electron droplet the curves show a dip, followed by a peak which 
does not change size or shape when the number of electrons changes, in contrast to the 
behaviour seen in Fig. |l|. 

Fig. ^ also demonstrates how the curves go to zero when the two quasi-holes overlap, 
i.e. r ^ 0. In Fig. |^ this small relative distance behaviour is shown on a larger scale for five 
different electron numbers, and compared to the Berry phases calculated from two different 
states in the system of two ideal anyons. The three lowest lying curves are, from bottom 
and up, for 75 electrons, then coinciding curves for 20, 50 and 100 electrons, and finally 
for 200 electrons. We conclude that the behaviour at small relative distance is due to the 
finite size of the quasi-holes themselves and is not affected by the number of electrons in 
the system. 

The two highest curves in Fig. |^ do not refer to the quantum Hall system, but offer an 
explicit comparison of features of the Laughlin quasi-holes to the case of ideal, point like. 



anyons, as discussed in the Appendix. In Ref. it was pointed out that the first order 
correction to the constant, asymptotic value of the two quasi-hole Berry phase favours the 
projected anyon position eigenstate rather than the coherent state of the su(l,l) algbra as 
the localized anyon state that best fits to the Laughlin quasi-hole state. This conclusion 
was based on the plasma analogy, saying essentially that the corrections to perfect screening 
should fall off exponentially rather than algebraically. It is confirmed by Fig. where the 
curve representing the coherent state lies highest, and the curve representing the projected 
anyon state indeed lies closer to the quasi-hole curves. 

The fact that the curves in Fig. ^ go to zero for r — > 0, can be understood in the 
anyon model as due to the impossibility of constructing states with the point like anyons 
completely localized, using only states belonging to the lowest Landau level. The parameter 
z, or z*, defining the anyon states we consider, is different from the relative position of the 
anyons, and moving z to —z around a circle of radius r is not the same as interchanging 
the anyon positions. Hence, in the limit r — 0, no anyon statistics phase is picked up. 

In Fig. 1^ we show the bulk behaviour of u^^, as calculated for 100 electrons. Again we 
take this as a confirmation, to at least three decimal digits, of the known result of 1/3 for 
the quasi-hole statistics parameter. 

Our final topic is what we have up to now referred to as the statistics parameter 
of the quasi-electrons, i.e. z/^*^ in Eq. (|I3|). Fig. ^ shows this quantity for the cases of 
20, 50, 75, 100 and 200 electrons. The 100 and 200 electron curves had to be cut at 
r = 8 and r = 6, respectively, because of numerical problems. This figure faces us with 
the question of whether the Laughlin quasi-electrons can be said to have a well defined 
statistics parameter at all. In our opinion, the term "statistics parameter" is justified only 
if the quantity denoted is independent of the distance between the two quasi-particles as 
long as they neither overlap nor are close to the edge of the electron droplet. We have seen 
that this criterion is well satisfied by the quantity u^^, Eq. (O), but Fig. S reveals that it 
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definitely does not hold for u'^'^, Eq. (|T3|), for the system sizes we have considered. 

We would like to stress that if we had neglected the extra terms in the integrals in 
Eqs. ( PBD and ( PBD as compared to Eqs. (^) and then this strange behaviour would 
not have appeared, and the statistics parameter for the quasi-electrons would have been 
equal to minus the statistics parameter of the quasi-holes. Thus, the fact that the computed 
z/'^'^ and have the same sign, is alone proof that the extra terms are essential. 

Fig. ^ shows the same curves on an enlarged scale. It also includes two curves that offer 
an explicit comparison to the case of ideal, point like anyons, for small r these lie below 
and to the right of the curves referring to the quantum Hall system. Note that the five 
quasi-electron curves depend only slightly on the number of electrons for small r, and thus 
offer a signature of the Laughlin quasi-electrons themselves. Note also that all the curves 
start out negative at small r before becoming positive. 

The two curves referring to the anyon system are calculated from states that are singular 
but normalizable. In fact, since the calculated right hand side of Eq. ( [T^ ) is positive, the 
only possible way to extract an anyon parameter is to make a correspondence with such 
singular states in the ideal anyon system. Fig. ^ shows that the coherent state is closest 
to the quasi-electron curves for small r, whereas the projected position eigenstate gives a 
better fit for large r. 

Fig. 1^ indicates that the bulk value of may possibly become well defined {i.e. 
constant) in the thermodynamic limit. But if so, it is not obvious that the limiting value 
would be 1/3. 



5 Concluding remarks 

To summarize, our numerical study has nicely confirmed the results earlier presented in 
the literature, that the quasi-holes at filling fraction 1/3 have charge e/3, and have anyon 
parameter and exclusion statistics parameter both equal to 1/3. 

If it were true that the normalization factors Ji and I2 were the same for the Laughlin 
quasi-electrons and quasi-holes, then the Berry phases for the quasi-electrons would be 
the same as for the quasi-holes, only with opposite signs. This would then lead us to the 
numerical results q^^ = —e/3 and u*^^ = —1/3, whereas the exclusion parameter would be 

uT = +1/3. 

However, for the quasi-electrons we have seen that for systems with up to 200 electrons 
the charge is affected by the extra terms discussed above, so as to give a bulk value that is 
slightly lower than the expected —e/3. The deviation is quite small and may well vanish 
for infinitely many electrons, but we do find it interesting to note that the quasi-electrons. 



which have a size of the order one magnetic length [|Tl|, seem to be affected by the edge 
even when it is as much as ten magnetic lengths away. 

For two quasi-electrons encircling one another our results differ dramatically from those 
earlier proposed in the literature. The extra terms affect the Berry phase such that it is 
unclear whether one can extract any statistics parameter at all, even in the thermodynamic 
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limit. What we can conclude is that the the statistics parameter, if at all meaningful, has 
a sign which is opposite to what it would have if the extra terms were unimportant. It 
may possibly converge to 1/3 in the thermodynamic limit, even though we find no strong 
indication that it does. If the Laughlin quasi-electrons have an anyon parameter u^'^ = 1/3, 
it means that the anyon wave functions are singular, and the exclusion statistics parameter 
is = —u^'^ = —1/3. This result fits in with the naive argument saying that if the quasi- 
holes and quasi-electrons are anti-particles, they should have the same anyon parameter. 
However, it is in not in agreement with the value 2 — 1/3 found from state counting, which 
corresponds to non-singular anyon wave functions. 
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A Berry phases of the two-anyon states 

In order to model the Laughlin quasi-particles as anyons of statistics parameter u in the 
lowest Landau level, it is necessary to establish a correspondence with specific states in 
the anyon system. There are two different two-anyon states that are natural to consider as 
states that maximally localize each of the particles |1T3| , p!8[] . They are the coherent state of 
the su(l,l) algebra and the anyon coordinate eigenstate projected onto the lowest Landau 
level. These states are identical in the fermion and boson cases, but not in general. 

Assume that the magnetic field is B and the anyon charge is e/m, with e > and 
B > 0. The relative angular momentum of the two anyons, in the lowest Landau level, 
have eigenvalues 2fc + z/ with k integer. Let \k, v) be the corresponding orthonormal eigen- 
states, and let \u) be the position eigenstates, with u the complex relative coordinate. 
Asymptotically as 1^1 — we have that 



Usually one requires that 2A; + > 0, so that the angular momentum eigenstates are non- 
singular. But if— l<2A; + z/<0 the corresponding state would be normalizable, although 
singular, and one may include it. 
The coherent state can be written 



{u I A;, z/) oc u 



2k+u 



(44) 



oo 



Z' 



.*2k 




(45) 



whereas the projected position eigenstate is 




(46) 
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Here A4 and M'^ are normalization factors, and 2; is a parameter labelling the states, such 
that 2z is to be interpreted as the relative position of the two anyons, measured in units 
of \/2 times the magnetic length l/y/eB. This scaling of distances corresponds to the 
convention used for the quantum Hall system in this paper. The only restriction on the 
statistics parameter v in these formulae is that z/ > 0, if the states are required to be 
non-singular, or that v > —1, if singular but normahzable states are accepted. But note 
that the sums leave out one or more of the non-singular angular momentum eigenstates 
when z/ > 2. 

For the statistics parameter as measured by the Berry phase we find the expression 



CAr) = -r'-r'^ln ''IZ (47) 




for the coherent state, with Iu-1/2 a modified Bessel function, and 

^ ^ m dr^ V«^r(2/ + z/+l); ^ ^ 

for the projected anyon position eigenstate. In both cases, the one-particle contribution 
— 2r^/m has been subtracted, and the asymptotic value as r — 00 is z/. 

These functions are plotted in Fig. H, for z/ = 1/m = 1/3, and compared to the statistics 
parameter found from the Laughlin quasi-hole states. 

The exclusion statistics parameter, which is given by the asymptotic behaviour of the 
two- anyon wave function as |m| — 0, is z/i = z/ for the coherent state as well as for the 
projected position eigenstate. 

In the case of anyons of negative charge — e/m, the relative angular momentum eigen- 
values are —2k + v with k integer, and the relative angular momentum eigenstates have 
the asymptotic form 

(u I fc, z/) oc (m*)2'="" (49) 

as \u\ — > 0. Thus, the state \k^v) is non-singular if 2A; — z/ > and normahzable if 
2k-iy> -1. 

The formulae (^) and (^) still hold with the substitutions p —v and z* — >■ z, and 
with the restriction on the statistics parameter v that z/ < 0, if only non-singular states 
are accepted, or that z/ < 1, if singular but normahzable states are accepted. A further 
consequence of the substitution z* z is a, change of signs in the Berry phase formulae 
(^) and (^81). The functions —C^y{r) and —A_y{r), with z/ = 1/3, are plotted in Fig. ^ 
for comparison to the statistics phase as computed for the Laughlin quasi-electrons. 

The exclusion statistics parameter, given by the asymptotic behaviour in the limit 
ImI 0, is now seen to be Pi = —p. 
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The quasi-hole charge (f^ je^ Eq. (|T0|), 
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as a function of r, the quasi-hole distance 
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Figure 3: The quasi-electron charge g^^'/e, Eq. (0), as a function of r, the quasi-electron 
distance from the origin. The figure presents results for, from left to right, 20, 50 and 75 
electrons. The horizontal line is —1/3. 
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Figure 4: Quasi-electron charge, (f^^ je. The number of electrons, from the lowest lying 
curve and up, is 20, 50, 75, 100 and 200. The highest curve is the constant —1/3. 
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Figure 5: The quasi-hole statistics parameter u"^^, Eq. (|T2D, as a function of r, half the 
distance between the two quasi-holes. The curves are, from left to right, for 20, 50, 75, 100 
and 200 electrons, and the horizontal line is 1/3. 
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Figure 6: Small relative distance behaviour of , and Berry phases calculated from two 
different localized states in the system of two anyons. The lowest lying curve is for 75 
electrons, then follows a common curve for the three cases of 20, 50 and 100 electrons, and 
the third curve is for 200 electrons. Somewhat higher lies the Berry phase curve calculated 
from the anyon position eigenstate projected onto the lowest Landau level, and even higher 
the one calculated from the coherent state of the su(l,l) algebra. See the Appendix for 
details. 
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Figure 7: Quasi-hole statistics parameter i/"^*^ for 100 electrons, compared to 1/3, empha- 
sizing the bulk behaviour. 
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Figure 8: The quantity u"^"^, Eq. (p!3|), to be interpreted as the quasi-electron statistics 
parameter, versus r, half the distance between the quasi-electrons. The lowest lying curve 
is for 20 electrons. Next are curves for 50 and 75 electrons. The 100 electron curve is cut 
at r = 8 and the 200 electron curve at r = 6, to avoid numerical problems. The horizontal 
line is 1/3. 
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Figure 9: Small relative distance behaviour of u'^^ for the cases 20, 50, 75, 100 and 200 
electrons. Also shown arc two curves corresponding to ideal anyons. The five curves for 
the different electron numbers coincide for small r. The 200 electron curve overshoots 
the horizontal line 1/3 at r fti 3. The curve lying lowest for small r represents an anyon 
eigenstate projected onto the lowest Landau level, whereas the curve going highest for large 
r represents the coherent state of the su(l,l) algebra. 



25 



